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A comprehensive theory is presented for the voltage, temperature, and spatial dependence of the 
tunneling current between a scanning tunneling microscope (STM) tip and a metallic surface with an 
individual magnetic adatom. Modeling the adatom by a nondegenerate Anderson impurity, a general 
expression is derived for a weak tunneling current in terms of the dressed impurity Green function, 
the impurity-free surface Green function, and the tunneling matrix elements. This generalizes 
Fano's analysis to the interacting case. The differential-conductance lineshapes seen in recent STM 
experiments with the tip directly over the magnetic adatom are reproduced within our model, as 
is the rapid decay, ~ 10A, of the low-bias structure as one moves the tip away from the adatom. 
With our simple model for the electronic structure of the surface, there is no dip in the differential 
conductance at approximately one lattice spacing from the magnetic adatom, but rather we see 
a resonant enhancement. The formalism for tunneling into small clusters of magnetic adatoms is 
developed. 



PACS numbers: 72.15.Qm, 61.16.Ch, 72.10.Fk 
I. INTRODUCTION 

Recently, a scanning tunnelingjjqicroscope (STM) was 
used in two separate experimentsEro to directly probe the 
local electronic structure of an isolated magnetic adatom 
on a metallic surface. By measuring the tunneling cur- 
rent from the substrate to the STM tip, a narrow res- 
onant feature was seen in the differential conductance 
when the tip was placed directly above the magnetic 
adatom: an antiresonance in the case of a cerium adatom 
on the (111) surface of silver ,EI and an asymmetric Fano 
resonance in-the case of a cobalt adatom on the (111) sur- 
face of gold.0 These resonant structures gradually disap- 
peared upon removing the tip from the adatom. Specif- 
ically, the asymmetric lineshape for Co on Au(lll) first 
evolved into a more symmetric dip, before disappearing 
altogether at a lateral distance of about 12 A from the Co 
adatom. □ Both experiments were interpreted as a mani- 
festation of the Kondo resonance that develops due to the 
screening of the local moment on the magnetic adatom 
by the substrate conduction electrons. 

While similar observations of the Kondo effect for a 
single magnetiCjimpurity were recently reported both in 
quantum dotsBu and in metallic point contactsjj STM 
spectroscopy offers the unique ability to spatially resolve 
the electronic structure around the magnetic adatom. 
This provides direct information about the screening of 
the local moment, allowing for critical comparison be- 
tween theory and experiment. Another novel aspect 
of STM spectroscopy in this context is the quantum- 
mechanical interference between direct tunneling into the 
magnetic adatom and tunneling into the underlying sutt 
strate conduction electrons, as first discussed by Fanccl 
for the noninteracting case. For a quantum dot placed in 



between two metallic leads, the analogous interference is 
between Kondo-assisted tunneling through the dot and 
direct tunneling between the leads, the latter being ex- 
tremely small. 

Motivated by the recent STM experiments, this paper 
provides a comprehensive theory for the voltage, temper- 
ature, and spatial dependence of the tunneling current 
between an STM tip and a metallic surface with an in- 
dividual magnetic adatom. Modeling the adatom by a 
nondegenerate Anderson impurity, a general expression 
is derived for a weak tunneling current in terms of the 
fully dressed impurity Green function and the impurity- 
free surface Green function. The impurity Green func- 
tion is evaluated in turn using the non-crossing approx- 
imation (NCA), while the surface Green function is ob- 
tained within a tight-binding model with free boundary 
condition at the surface. This allows for a consistent de- 
scription of the energy and the spatial dependence of the 
tunneling current, as is required in this case. Both cases 
of point tunneling and that of a finite spatial extent in 
the tunneling matrix element between the tip and the 
substrate conduction electrons are considered within this 
model. 

As expected of the Kondo effect, a sharp resonant 
structure is found to develop in the low-voltage differ- 
ential conductance as the temperature is lowered down 
to the Kondo temperature, Tk- The width of the reso- 
nance is proportional to Tk at low temperatures, grow- 
ing with T for T > Tk ■ Its shape is governed by a single 
interference parameter a, much in the same way as in 
the noninteracting caseB The value of q and its spatial 
variation depends quite sensitively on details of the un- 
derlying surface Green function and the tunneling matrix 
elements, which limits the direct applicability of our re- 
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suits to the experimental data. Nevertheless, we are able 
to make some qualitative statements regarding the ex- 
periment. In particular, the indirect interference with 
the magnetic adatom, i.e., that due to the tunneling be- 
tween the STM tip and the underlying conduction elec- 
trons, is found to have a characteristic range of the order 
of two lattice spacings, consistent with the limited spa- 
tial extent of the Kondo resonance in the experiment. In 
addition, the spatial dependence of the differential con- 
ductance as measured for Co on Au(lll)a is shown to be 
nongeneric, indicating a crucial dependence on the micro- 
scopic details of the (111) surface of gold. Finally, since 
the widths of the resonant features in Refs. [I] and |2| are 
both notably larger than the temperature, we conclude 
that T < Tk in both experiments. 

The remainder of the paper is organized as follows: In 
Sec. H the basic model is introduced, and a general ex- 
pression is derived for the differential conductance in the 
limit of a weak tunneling current. Sec. Ill then details our 
calculation scheme for the dressed impurity Green func- 
tion and the impurity-free surface Green function. Sim- 
plified expressions for the first few nonlocal surface Green 
functions of our model are provided in Appendix |A|. Re- 
sults for the voltage, temperature, and spatial depen- 
dence of the low- voltage differential conductance are pre- 
sented in turn in sections |v| and with Sec. |v| ded- 
icated to the case of point tunneling between the STM 
tip and the underlying substrate conduction electrons, 
and Sec. |v| devoted to the case of a spatially extended 
matrix element for tunneling. The main results of the 
paper are finally summarized in Sec. 



H tip =)J. S ub + eV 



VI 



followed by a 

brief comparison with previous work in Appendix and 
a generalization of our approach to the case of multiple 
magnetic adatoms in Appendix |c]. 



II. MODEL AND TUNNELING CURRENT 

The system under consideration is shown schematically 
in Fig. [IJ. It consists of an individual magnetic adatom, 
modeled by a d level with energy ed and an on-site re- 
pulsion U, deposited on top of a metallic surface. The 
d electrons hybridize locally with the underlying conduc- 
tion electrons via the matrix element Vh- This setting 
is probed by an STM tip which is placed directly above 
the surface point R s , and which couples separately to the 
d electrons and to the local conduction electrons at site 
R s through the tunneling matrix elements td and t c , re- 
spectively. The adatom is taken to be positioned at Ri. 
An applied voltage bias between the substrate and the 
tip offsets the two chemical potentials, [itip — Hsub = eV, 
causing a weak electrical current to flow between the sub- 
strate and the tip. 

Setting the substrate chemical potential as our refer- 
ence energy, the corresponding Hamiltonian takes the 
form H = Tisub + Hup + Htun, where 




FIG. 1. Schematic description of the physical system. A 
metallic surface with an individual magnetic adatom is probed 
by an STM tip. The tip couples separately to the atomic d 
electrons and to the underlying conduction electrons via the 
tunneling matrix elements td and t c , respectively. A voltage 
bias is applied between the sample and the STM tip, causing 
a weak electrical current to flow between the substrate and 
the tip. 
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Here ct (at ) creates a bulk (tip) conduction electron 

kcr V kcr' K 1 ' 

with wave number k and spin projection a, and d\ creates 
an atomic d electron with spin a. The fermion operators 
ipcr{r) and A a in Eqs. (@)-(§) are the local conduction 
electrons at point r on the surface and at the edge of the 
STM tip, respectively. Explicitly, 



k 



(4) 
(5) 



where tp%(r) an d are the corresponding bulk and tip 
single-particle wave functions, evaluated at point r on 
the surface and at the edge of the tip, respectively. The 
atomic energies U+e d and — ed are assumed to be positive 
and large, such that a local moment forms on the adatom. 

For zero tunneling, t c ,td — 0, the local d moment 
undergoes Kondo screening below a characteristic tem- 
perature Tk oc exp [— 1/poJ], where po is the local 
surface density of states at the Fermi level, and J = 
2V£ [l/\ed\ + 1/(U + eo)] is the effective exchange cou- 
pling between the d moment and the underlying ct con- 

_ kcr 

duction electrons. For nonzero tunneling and a finite 



2 



voltage bias, ane is dealing in principle with a nonequi- 
librium state. Q'u However for a weak tunneling current, 
as is the case in the experiments of Ref. [j] and Ref. [|, 
the differential conductance essentially probes the zero- 
tunneling (equilibrium) Kondo resonance. To see this we 
note that, for weak tunneling, it is sufficient to evaluate 
the tunneling current from the substrate to the STM tip 
to second order in t c and t^. Using standard diagram- 
matic techniques one obtains 

rj /-OO 

1{ y) = -* E / p^ e - eV )pf^ - eV )-fw 

™„ J-oo 

(6) 

where /(e) is the Fermi-Dirac distribution function, and 
PAa{t) = -ImG^„(e + irj) and pfa(e) = -ImG/ (7 (e + if]) 
arc the zero-tunneling spectral functions corresponding 
to A a and 



fa = t c 1p,j(R s ) + t d d a 



(7) 



respectively. Here we have introduced the notation by 
which the arguments of G Aa and Gf a (and thus those 
of pAa and pfa) are measured relative to their respective 
chemical potentials, i.e., p,u p — eV and p, S ub = 0. 

All information of the Kondo effect in Eq. (^|) is con- 
tained within the retarded Gf a propagator, 

Gfa (e + ill) = t 2 c G &3 A (e + irj) + G d a (e + irj) x 
td + tcVhG^^e + irj) t d + t c V h G^. n3 (e + ir]) , (8) 

which features the fully dressed impurity (d-electron) 
Green function in the absence of tunneling, G^(e + irj), 
and the impurity-free surface Green function, 



Gf u f 2 (e + ir]) = ^ 



f^i)^*^) 
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The differential conductance is given in turn by the 
derivative of Eq. (^) with respect to V. Assuming 
PAais) — Pa is essentially energy independent on the 
scale of the voltage bias and the temperature, one ob- 
tains 



G{V) 



2f_ 



Pa 



:p f a(e + eV)^-de. (fO) 



EquationsJ|8])-(|Trj) are a generalization of Fano's orig- 
inal analysiso to the interacting case. They express the 
differential conductance in terms of the fully dressed im- 
purity Green function, the impurity-free surface Green 
function, and the microscopic tunneling parameters of 
the model. The remainder of the paper is devoted to 
evaluation and analysis of these equations, starting with 
the Green functions Gf(e + irj) and Gp ly p 2 (e + irj). 



III. IMPURITY AND SURFACE GREEN 
FUNCTIONS 

A. Impurity Green function 

We have evaluated the impurity Green function in the 
limit U — * oo, where double occupancy is forbidden on 
the magnetic adatom. Calculations were performed using 
the non-crossing approximations (NCA) , which is a self- 
consistent perturbation theory about the atomic limit. 
This approach is known to provide a good quantitative 
description of the temperature range T > Tk, and has 
the advantage of being capable of treating systems with 
realistically small Kondo temperatures. This is to be con- 
trasted with Quantum Monte Carlo approaches, which 
can not go to very low temperatures. As a large- N the- 
ory, however, the NCA fails_to recover Fermi-liquid be- 
havior at low temperatures and it overshoots the uni- 
tary limit for T < Tk in the N — 2, nondegenerate case. 
To avoid these pathologies of the NCA, we restrict atten- 
tion in this paper to the range T > Tk- 



B. Surface Green function 

Due to the local nature of the hybridization between 
the d and conduction electrons in Eq. (|l|), only the local 
conduction-electron density of states is relevant to the 
impurity Green function, G^(e+irj). Conventionally, this 
allows one to parameterize the band by a single function 
for the density of states, often chosen for convenience to 
have a Lorentzian, box, or Gaussian form. In contrast, 
calculation of the differential conductance for R s ^ Ri 
requires detailed information of the band dispersion. In 
particular, a consistent theory for the energy and spatial 
dependence of G(V) requires one to start from a micro- 
scopic description of the underlying conduction band. 

As a generic model for the substrate conduction band, 
we consider a simple-cubic tight-binding Hamiltonian, 
with open boundary conditions at the surface. Setting 
the lattice spacing as our unit length, the lattice is de- 
scribed by the integer grid r*j = (xi,yi,Zi) with Zj > 0, 
such that the first (surface) monolayer corresponds to 
z.i = 1. In each of the x and y direction we take peri- 
odic boundary conditions, i.e., Xi and yi are equivalent 
to Xi + N and + N, respectively. In the z direction 
we apply open boundary conditions, namely, cp ia (the 
conduction-electron annihilation operator at site r*) van- 
ishes for Zi — and Zi = N. The corresponding Hamil- 
tonian reads 



T~Cbar. 



t 



E 

<i,j>,a 



fct 



•ct 



(ii) 



where <i,j> denotes nearest-neighbor lattice sites. For 
Psub = 0, the case considered here, the Hamiltonian of 
Eq. (O) describes a half-filled band. 
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The tight-binding Hamiltonian of Eq. ( |ll| ) is diagonal- 
ized by converting to the single-particle basis 



AT3 X/ 



-i(k x Xi-\-k y yi) 



&hi(k z Zi 



(12) 



Here k takes the values (2n x ,2n y , n z )ir/N, with the inte- 
gers < n x , n y < N and 1 < n z < N. The corresponding 
single-particle energies take the form 



££ = 2t ^ cos(fci). 

l=x,y,z 



(13) 



Since the tight-binding wave functions ip%{f) are con- 
fined to a discrete set of lattice sites, the Green function 
of Eq. (j^) is meaningful only in the case where both r\ 
and f*2 are lattice sites. Taking the thermodynamic limit, 
N — > oo, one thus obtains 

2 r r r 

Gfi,P 2 ( e + W = —3 / dk x / dk y dk z (14) 
n Jo Jo Jo 

sin(fc z zi) sm(k z z 2 ) 



x cos[k x (xi-x 2 )] cos[ky(yi-y 2 )}- 



e k + l7 l 



from which the surface Green function follows by setting 

z\ = z 2 = 1: 

2 

Gi m (e + irj) = —s I dk x / dk y / dk z 
7T Jo Jo Jo 

s'm 2 (k z ) 



x cos(k x l) cos(k y m)- 



e fc + "? 



(15) 



Here and in the remainder of the paper we label the sur- 
face Green function by two integers: I — x\ — x 2 and 
m — Vi ~ V2- We further note that Gi m (e + irf) and 
G m i (e + irf) are identical by symmetry, and are only de- 
pendent upon the distances \l\ and \m\. 

Equation ([l^) can be further simplified to just a single 
integration by exploiting the relation to the local tight- 
binding Green, function for a two-dimensional, simple- 
cubic lattice £3 



G 2D (0 - -K[(4t/Cn 



(16) 



Here K(() is the complete elliptic integral of the first 
kind, axialytically continued to the upper and lower half 
planes.Eil Specifically, the local surface Green function, 
Goo, is conveniently expressed as 

2 C 

Goo(e + iff) = — / dk sin 2 (k)G 7 d[c — 2t cos(/c) + ir]], 
n Jo 

(17) 

with analogous expressions for the first few nonlocal 
Green functions (see Appendix |A|) . 

Figure ^displays the first four surface Green functions. 
Both the real and imaginary parts of G; m (e + irf) decay 
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FIG. 2. Real and imaginary parts of the tight-binding sur- 
face Green function G; m (e + irf), for different values of (l,m). 



as a function of y/Y 2 + m 2 , alternating between even and 
odd functions of energy. For even (odd) |Z| + \m\, the 
real part of Gi m {e + irf) is odd (even) in energy, while the 
imaginary part is even (odd). As a result, G; m (0 + irf) 
is purely imaginary when \l\ + \m\ is even, and purely 
real when \l\ + |m| is odd. This has a profound effect 
on the spatial dependence of the low-voltage differential 
conductance for this model, as a purely real G; m (0 + 
irj) gives rise to resonant enhancement of the differential 
conductance for R s — Ri = (l,m,0). By contrast, an 
imaginary Gi m (0+irj) can result both in an antiresonance 
or an asymmetric Fano resonance, depending on details 
of the tunneling matrix elements t c and td- We shall 
return to this point in great detail below. Finally, we 
note that the imaginary part of Gi m (0 + irf) is nonzero 
only in the range —6t < e < 6t, and that the van-Hove 
singularities at e = ±2t are considerably smoother at the 
surface than deep inside the bulk. 



IV. RESULTS FOR POINT TUNNELING 

We have evaluated the differential conductance accord- 
ing to Eqs. (P)-(]lO|), within the tight-binding model of 
Eq. ( pd| ) for the underlying conduction band. if^Ri) 
was identified in this picture with the creation of an elec- 
tron at the lattice site Ri = (0,0,1). The local density 
of states used in computing the impurity Green function 
was taken accordingly to be p{e) = — ilm {Goo(e + 
where Goo is given by Eq. (|17|). The corresponding den- 
sity of states at the Fermi level is equal to po = p(0) = 
0.525/t. Focusing on the case where the STM tip is po- 
sitioned directly above the lattice site R s — (l,m,l), 
each of the surface Green functions G^ 3 (e + irf) and 
Gr. r ( e + * ? 7) m Eq. (Sh was identified with Gi m {e + irf). 
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FIG. 3. The differential conductance, G(V), as a function 
of voltage, for an STM tip placed directly above the mag- 
netic adatom. Here T = Tk, ta = 0, and Go = 4e 2 t^.popA/h. 
The impurity model parameters are Ed/t = —1.67, T/t = 0.4, 
and U — oo, corresponding to a Kondo temperature of 
Tk /t — 10 . The impurity Green function, whose spectral 
part Af(e) = — (l/7r)Im {G^(e + z7/)} is plotted in the inset, 
was computed using the NCA. In addition to the standard 
broad feature near the d level, a sharp antiresonance is seen 
in G(V) at zero bias, corresponding to the Abrikosov-Suhl 
resonance in the impurity spectral function. 



A. Tip placed above the magnetic adatom 

We begin with an STM tip positioned directly above 
the magnetic adatom. As the temperature is lowered 
down to Tff, a sharp resonant structure develops in the 
differential conductance at zero bias, in addition to the 
standard broad feature near the d level. This is demon- 
strated Fig. | for T = T K and t d = 0. The sharp an- 
tiresonance seen at zero bias in this case directly corre- 
sponds to the Abrikosov-Suhl resonance in the impurity 
spectral function, which is plotted for comparison in the 
inset of Fig. |^. Here and throughout the paper, the im- 
purity Green function was calculated within the NCA, 
using the impurity model parameters e^/t = —1.67, 
T/t = npoVf/t = 0.4, and U = ooH The correspond- 
ing Kondo temperature, Tjc/t = 10™ 3 , was defined as 
the temperature at which the "resistivity" integralJ13 



R(T) = 



dm 



lm{Gi(e + ir])} de 



de 



(18) 



reduces to 50% of its zero-temperature limit. This def- 
inition of Tk agrees to within a factor of order unity 
with the half-width of the Abrikosov-Suhl resonance in 
the impurity spectral function. 

Focusing hereafter on the low- voltage resonant struc- 
ture, Fig. depicts the resonance's dependence on go = 
td/ftpotcVh, which plays the role of Fano's interference 
parameter q in this case. Here G = Ae 2 t^p PA/^ is the 
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FIG. 4. The low-voltage differential conductance for an 
STM tip placed directly above the magnetic adatom. Here 
T = Tk and go = td/irpotcVh, with all other model parame- 
ters as in Fig. |j| With increasing |go|, G(V) evolves from an 
antiresonance for |go| "C 1, to an asymmetric Fano resonance 
for intermediate |qo|, to a full resonance for |go| 3> 1 (not 
shown) . 



zero-temperature conductance in the absence of an im- 
purity. With increasing \qo\, the low- voltage differential 
conductance evolves in Fig. ^ from an antiresonance for 
| go I ^ 1 , to an asymmetric Fano resonance for intermedi- 
ate | g 1, to a full resonance for |go| 3* 1, much in the same 
way as in the noninteracting case. The antiresonance for 
qo = resembles that for a Ce adatom on Ag(lll),El sup- 
porting the interpretation of Li et al. that the observed 
lineshape stems from Kondo screening of the Ce moment 
with only weak direct coupling between the tip and the 
adatom. The go = 1 curve is similar in turn to the Fano 
resonance observed by Madhavan et al. for a Co adatom 
on Au(lll),0 while larger values of go display too shal- 
low a dip and too high a peak as compared to the data 
of Ref. |. 

Figure || shows the temperature dependence of G(V), 
for the two representative values of go = and go = 0.7. 
The main effect of a temperature is to broaden and smear 
the resonant structure in G(V), whose width grows ac- 
cording to T for T > Tk ■ Notice the significant reduction 
in the overall resonance height by the time T/Tk ~ 10. 
This strong temperature dependence of the differential 
conductance has two contributions: (i) A rapid decrease 
in the Abrikosov-Suhl resonance in the impurity spectral 
function with increasing temperature, and (ii) Further 
smearing of the impurity Green function due to the con- 
volution with the derivative of the Fermi-Dirac function 
in Eq. (|l(]). The fact that the well-developed features in 
Refs. [I] and have widths that are notably larger than 
T is thus a clear indication that T < Tk in both experi- 
ments. 
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FIG. 5. Temperature dependence of G(V), for an STM tip 
placed directly above the magnetic adatom. Here qo = and 
qo — 0.7 in the lower and upper graphs, respectively. All other 
model parameters as in Fig. |§| 



B. Spatial dependence of the differential 
conductance 

As stated above, the novel aspect of STM spectroscopy 
of the Kondo effect is the ability to spatially resolve the 
electronic structure around the magnetic adatom. In 
Fig. ^ we have plotted the low-voltage differential con- 
ductance as a function of (l,m), for the idealized case 
where the STM tip is placed directly above the lattice 
site R s = (1,171, 1). Here each set of curves corresponds 
to a different fixed value of go = td/Trpot c Vh, which no 
longer corresponds for R s ^ Ri to Fano's interference 
parameter q. The latter also depends on the real and 
imaginary parts of G^ g (0 + irf), as discussed below. 
Physically one should recall, though, that qo oc de- 
cays to zero with increasing distance from the magnetic 
adatom, hence G(V) always tends to the go = curve as 
the lateral distance from the adatom is increased. 

Fixing the value of qo for the time being and increasing 
the lateral distance from the magnetic adatom, the low- 
bias differential conductance approaches the asymptotic 
curve 



G(V) = Go 



l + r q lJ°\m{Gi(e + ir,)}^-de 



(19) 



This result stems from the decaying nature of Gi m (0+ir]), 
along with the fact that Goo(e + if}) = —inpo is approxi- 
mately constant on the scale of the voltage bias and the 
temperature. As seen in Fig. [| Eq. (|l9| ) is approached 
at a lateral distance of about two lattice spacings from 
the magnetic adatom. For larger distances the resonance 
height is basically proportional to the residual coupling 
to the magnetic adatom squared. Specifically, there are 
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FIG. 6. Spatial dependence of the low-voltage differen- 
tial conductance, G(V), for T — Tk and different values of 
qo = td/irpotcVh- All impurity model parameters are as in 
Fig. |§| The individual curves are offset by one unit each, 
according to the lateral distance from the impurity site. 



no visible traces of the Kondo resonance for qo = when 
the tip and the adatom are two lattice spacings apart [the 
curve (I, m) — (0, 2) in Fig. |(|. Such a limited spatial ex- 
tent of the Kondo resonance in G(V) is consistent with 
the one seen experimentally for Co on Au(lll),B indicat- 
ing that qo is effectively zero in the experiment above a 
lateral distance of about 10A. Indeed, the low-bias reso- 
nant structure is also mostly gone by lateral distance of 
lOAfor Ce on Ag(lll)Jd even though the Ag(lll) surface 
state at — 70meV does not fully set in before a distance 
of about 40 A. 

The resonant enhancement of G(V) at a lateral dis- 
tance of one lattice spacing, which occurs for any value 
of go j is not seen in the experiment. This feature of the 
calculated differential conductance is traced back to the 
fact that Goi(0 + irj) is purely real for our tight-binding 
model, resulting in a differential conductance G(V) that 
is once again given by Eq. (fly), but with 



go -> go + 



-Re {G O i(0 +irf)} = g - 0.372 (20) 



[see Eqs. (||) and ([To|)]. Thus, irrespective of the actual 
go that applies to a lateral distance of one lattice spacing 
in the experiment of Madhavan et al., our present model 
fails to recover the dip-like structure seen experimentally 
at such a distance for Co on Au(lll).cl 

That Goi (0 + irj) is purely real is a generic feature of 
half-filled nearest-neighbor tight-binding models on bi- 
partite lattices. It is lost, however, when the system is 
away from half filling, or upon inclusion of a next-nearest- 
neighbor hopping term. This suggests that the dip-like 
structure seen experimentally up to a distance of ~10 A 
from the Co adatomEl is due to a qualitative difference 
in the underlying surface Green function for Au(lll). 
Indeed, the (111) surface of gold is known to have a her- 
ringbone reconstruction with regions of fee and hep or- 
dering,li-3 which differs from the simple-cubic structure 



G 



considered here. This may support the interpretation of 
a qualitative difference in the underlying surface Green 
function. 

It should be noted, though, that such an interpretation 
for Au(lll) demands that Im{G^ g (0 + ir))} is never 
small compared to irpoqo + Re{G^ a (0 + i?y)} in the rel- 
evant range in R s — R\, In particular, this rules out any 
oscillatory behavior of Im{G^ ^ (0 + ir/)} as a function 

of R s — Ri in this range. A brief examination of Fig. || 
reveals that this condition is actually quite restrictive, at 
least for the tight-binding model of Eq. ([ll]): Only in a 
narrow window of band fillings is Im{G; m (0 + ir])} con- 
sistently negative for all |/| -I- |m| < 2. We thus conclude 
that the spatial dependence measured by Madhavan et 
al. is certainly not generic, but depends on details of the 
underlying band. 



V. FINITE SPATIAL EXTENT OF THE 
TUNNELING MATRIX ELEMENT 

Thus far we have considered an idealized point tunnel- 
ing between the STM tip and the substrate conduction 
electrons at point R s . In practice, however, the tun- 
neling matrix element has a finite spatial extent about 
R s , which is reflected in the different lineshapes that are 
observed when the tijxis removed from the Co adatom 
in opposite directions.El On the level of the model, a spa- 
tially extended tunneling matrix element is accounted for 
by replacing ip a (R s ) in the first term of Eq. (||) with a 
weighed sum over the conduction-electron degrees of free- 
dom around R s : 
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FIG. 7. The differential conductance, G(V)/Go, as a func- 
tion of wi/wq, for T = Tk, qo = 0, and an STM tip 
placed directly above the magnetic adatom. Here both wo 
and wi are assumed to be real and positive. Note that Go 
itself varies as a function of wi/wo [see Eq. ([m])], taking 
the values Goh/4e 2 t 2 c p p A = 1,0.838,0.598, and 0.356 for 
wi/wo = 0, 0.25, 0.5, and 1, respectively. All impurity model 
parameters are as in Fig. |j| As the ratio Wi/wo is increased 
from zero to one, G(V) evolves from an antiresonance to a 
resonance, corresponding to an increase in the effective inter- 
ference parameter q [see Eq. fl27|) and accompanying text]. 



Accordingly, the zero-temperature conductance in the 
absence of an adatom is equal to 



G = - 



4e 2 tj 



/ wpwp, G^ 



Wf4>cr(Rs + f). 



Here we use the convention 



= 1, 



(21) 



(22) 



which fixes the separation of the local tunneling matrix 
element at each r into tp = t c wp. For the tight-binding 
model of Eq. ( pi] ) , the sum over r in Eq. ( pi| ) extends 
over those lattice sites close to the STM tip. 

Upon substituting Eq. (21) into the tunneling Hamil- 
tonian of Eq. (|3|), the differential conductance remains 
given by Eqs. (]q) — (|To|) , but with the following modifica- 
tions to Eq. (|): 



G R 3 A ^ W ? w ?> G R,+F.g..+i>" 
GRsA ™* a2 WpG R.+r.R\ ' 

T 

G Ri,R s ~~ * ^ W r G Ri.R.+r- 



(23a) 
(23b) 
(23c) 



(24) 

which properly reduces to Go = 4e 2 ^poP J 4/^ in the limit 
of point tunneling. 

To examine the effect of a finite spatial extent in the 
tunneling matrix element, we go back to the tight-binding 
model of Eq. (|ll| ) , and to the case where the tip is placed 
directly above the magnetic adatom, i.e., R s — (0,0,1). 
In addition to tunneling between the tip and the lattice 
point R s , we introduce a nonzero tunneling matrix ele- 
ment to each of the four surface nearest-neighbors of R s . 
Restricting attention to the isotropic case, one is left with 
two different wp parameters: wq for the tunneling matrix 
element to the lattice point R s , and wi for the tunneling 
matrix element to each of its four surface nearest neigh- 
bors. The normalization condition, Eq. (E27), then reads 



(25) 



Figures [?] and || depict the evolution of the low- 
voltage differential conductance G(V)/Gq as a function 
of wi/wq, for the two representative values of qo = and 
<7o = 0.7. Here we have focused for simplicity on the 



7 



CD 
> 
CD 



1.5 



1.0 



0.5 



0.0 



1 1 




q = 0.7 






X ' 
V' 

1 


If-''' 

\i 

1 

w,=0 

w =4w, 

w D =w, 



-0.05 



0.00 

eV/t 



0.05 



FIG. 8. Same as Fig. 0, but for go = 0.7. As the ratio 
wi/wo is increased from zero to one, G(V) evolves from an 
asymmetric Fano resonance to an antiresonance, correspond- 
ing to a decrease in the effective interference parameter q [see 
Eq. ( [27] ) and accompanying text]. 



case where wq and w± are both real and positive, yet the 
qualitative picture does not depend on this choice. The 
dramatic effect that a nonzero w\ has on the differential 
conductance in this case can be understood within Fano's 
interference picture. Substituting Eqs. ( |23| ) into Eq. (||) 
and taking wq and w\ to be real, the effective interference 
parameter a la Fano is given for R s — Ri by minus the 
ratio of the real and imaginary parts of 

t d + t c V h [w Q G O0 (0 + it]) + 4 Wl G i(0 + irj)] . (26) 

Here Eq. ( p6|) corresponds to the resulting expression in 
each of the square brackets of Eq. (@). Using Goo(0 + 
irj) = —itt pa and Goi(0 + if]) = — 0.3727r j oo, as is appro- 
priate for the model of Eq. ( pi] ) , one obtains 



q=(q - 1A9wi)/w . 



(27) 



For g = 0, Eq. ( gTj ) reduces to —lA9(w 1 /w ), which 
varies from q = to q — —1.49 in going from w% = to 
wi = wq. This strong change in q produces the transition 
from an antiresonance to a resonance in the differential 
conductance of Fig. [?]. Similarly for qo = 0.7, Eq. (|2^ ) 
varies from q = 0.7 to q = 0.075 in going from wo = 1 to 
wo = w\ , which causes the transition from an asymmetric 
Fano resonance to an antiresonance in Fig. |^. 

Repeating the same argumentation for the case where 
R s and Ri are one lattice spacing apart, i.e., R s — 
(±1,0, 1) or R s = (0,±1, 1), Eq. @ is modified to 



td + t c Vh [w Gqi + wiGqq + 2wiGn + wiG 



02 J 



(28) 



with all Green functions evaluated at zero frequency. Us- 
ing Gn (0 + irj) = 0.354i7rp and G 02 (0 + irf) = 0.097iwp 
one obtains 
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FIG. 9. The differential conductance, G(V)/Go, as a func- 
tion of wi/wo, for an STM tip one lattice spacing removed 
from the magnetic adatom [i.e., R s — (±1, 0, 1) or (0, ±1, 1)]. 
As in Fig. |^, wo and u>i are assumed to be real and positive, 
with all impurity model parameters the same as in Fig. ^. 
For both qo = and go = 0.7, the differential conductance 
continues to feature a resonance for all < w\/wq < 1. 



which is large in magnitude throughout the range < 
w\ < wq, for both go = and go = 0.7. Given the large 
value of |g|, one expects the differential conductance to 
continue to show a resonance for both values of go and all 
< Wi/ Wo < 1, which is precisely what is seen in Fig. ||. 

Thus, while the inclusion of nonlocal tunneling be- 
tween the STM tip and the underlying substrate con- 
duction electrons obviously increases the parametric de- 
pendence of g, it does not necessarily assist in producing 
a dip in the low-voltage differential conductance when 
the tip and the adatom are one lattice spacing apart. 



VI. DISCUSSION 

We begin our discussion with the case of an STM tip 
placed directly above the magnetic adatom, depicted in 
Figs. ||-[| and Figs. fj]-|^. Similar to the noninteract- 
ing case, the shape of the Kondo resonance in the low- 
temperature, low-voltage differential conductance is gov- 
erned by a single interference parameter g, which depends 
both on the ratios of the tunneling matrix elements and 
on details of the impurity-free surface Green function. 
Specifically, in Fig. | q is equal to go = td/npotcVh, 
whereas in Figs. 0-||it is modified according to Eq. fl27|). 
The effect of a temperature is to rapidly broaden and 
smear the Kondo resonance in G(V), whose width grows 
according to T for T > Tk- This behavior stems both 
from the standard convolution with the derivative of the 
Fermi-Dirac distribution function in Eq. ([To]), and from 
the rapid decrease in the Abrikosov-Suhl resonance with 



8 



increasing temperature. Indeed, the qualitative differ- 
ence in the low-temperature, low-voltage differential con- 
ductance for a magnetic adatom as compared to that of 
a conventional noninteracting resonance is in the strong 
energy and temperature dependence of the impurity self- 
energy, whose real and imaginary parts cannot be re- 
garded a constant. 

Expecimentally, the antiresonance observed for Ce on 
Ag(lll)lil is similar to the go = curve in Fig. ||. As- 
suming point tunneling between the tip and the un- 
derlying conduction electrons, this supports the inter- 
pretation of weak direct tunneling between the tip and 
the Ce adatomji which requires at the same time that 
Gfi * (0 + ir)) is mostly imaginary. While plausible, this 
scenario is certainly not exclusive, as there are various 
other ways to obtain q « without resorting to a neg- 
ligible coupling to the adatom. For example, the curve 
wo = wi in Fig. || (corresponding to a spatially extended 
tunneling) is also characterized by |g| <C 1, although 
go = 0.7 is by no means small. Likewise.while the sim- 
ilarity between the g = 1 curve in Fig. and the Fano 
resonance for Co on Au(lll)0 is suggestive of compara- 
ble contributions from the tunneling to the adatom and 
to the underlying conduction electrons, one cannot rule 
out other combinations for which qo ~ 0. For example, 
setting qo ~ and wq/wx w —1.5 in Eq. (^7j ) also results 
in q w 1. On the other hand, the fact that the well- 
developed features in Refs. [l] and || have characteristic 
widths that are considerably larger than the temperature 
is a clear indication that T < Tr in both experiments. 

More detailed information about the underlying elec- 
tronic structure is contained in the spatial variation of 
the differential conductance, as measured, for example, 
in Ref. || for a Co adatom on Au(lll). Here, although we 
have considered a particular tight-binding model for the 
impurity-free surface, there are some qualitative state- 
ments we can make with regard to the experiment. Pri- 
marily, as seen in Fig. the indirect interference with 
the magnetic adatom, i.e., that due to the tunneling be- 
tween the STM tip and the underlying conduction elec- 
trons, is suppressed above a lateral distance of about two 
lattice spacings from the adatom. Hence the characteris- 
tic range for the indirect interference with the adatom is 
of the order of two lattice spacings. While this range may 
certainly depend both on details of the underlying band 
and on the presence of nonlocal tunneling between the 
tip and the substrate conduction electrons, we expect a 
qualitatively similar result for other microscopic models. 
From the limited saatial extent of the Kondo resonance 
for Ce on Ag(lll)U and Co on Au(lll)H we thus con- 
clude that go is effectively zero above a lateral distance 
of about 10A in these experiments. 

Our calculations further indicate that the spatial de- 
pendence of the differential conductance as measured for 
Co on Au(lll) is not generic, but intimately depends on 
the microscopic details of Au(lll). Indeed, while Madha- 
van et al. observe a dip-like structure that persists up to 
a lateral distance of ~10A from the adatomja we typically 



find a resonance at a distance of one lattice spacing. In 
Fig. [(], where point tunneling is assumed, this resonance 
occurs for any value of go, which is a special feature of 
the half-filled nearest-neighbor tight-binding model used. 
One may anticipate, though, a similar resonance within 
a range of one lattice spacing for other lattice models at 
half filling, since « (0 + iff) is expected to oscillate 
on a length scale of one lattice spacing. We emphasize, 
however, that the spatial dependence measured by Mad- 
havan et al. remains quite restrictive for the Hamiltonian 
of Eq. (|ll|) both away from half filling and in the case 
of nonlocal tunneling between the tip and the substrate 
conduction electrons. 

The particular tight-binding model used in this paper 
clearly limits the application of our results to the exper- 
iments. To make direct contact with the experimental 
data it is necessary to employ realistic Green functions 
for the (111) surfaces of silver and gold, which may be ob- 
tained, for example, from ab-inito calculations. It would 
be interesting to see if the combination of ab-inito calcu- 
lations for the (111) surface of gold with NCA calcula- 
tions for the many-body Kondo resonance can reproduce 
the particular spatial dependence of the differential con- 
ductance as seen for Co on Au(lll). 

Another interesting issue is the magnetic-field de- 
pendence of the low-temperature, low-voltage differen- 
tial conductance. With increasing magnetic field, the 
Abrikosov-Suhl resonance is first split for a moderate 
magnetic field, H ~ Tk, before a large magnetic field 
suppresses the Kondo effect altogether. A similar pat- 
tern is expected for the Kondo resonance in the differ- 
ential conductance. Unfortunately, treatment of a finite 
magnetic field within the NCA is hampered by the NCA 
pathology,EI hence a different approach is required. One 
possibility might be Quantum Monte Carlo simulations 
in combination with the Maximum Entropy method for 
analytic continuation,!^ although such an approach is 
restricted in treating realistically small Kondo temper- 
atures. 

Finally, in this paper we have focused on the case of an 
individual magnetic adatom; however, using the STM tip 
to atomically manipulate individual adatoms into form- 
ing small clusters, it might be possible to address the 
subtle interplay between the Kondo effect and magnetic 
correlations among the different adatoms. Most notably, 
the competition between the Kondo effect and antiferrtt. 
magnetic locking in the case of two close-by adatoms.EZl 
A first study of a Co dimer on Au(lll) along these lines 
was recently reported in Ref. |l8|. As detailed in Ap- 
pendix our formulation of the tunneling current is nat- 
urally extended to the case of multiple magnetic adatoms. 
Specifically, the single d-electron Green function enter- 
ing Eq. (@) is replaced by a matrix propagator, corre- 
sponding to all possible propagations within the adatom 
cluster. In this manner, one can analyze complicated 
multiple-adatom configurations in terms of the intra-site 
and inter-site d Green functions. 







It is our hope that the approach developed in this pa- 
per will prove useful in analyzing future STM measure- 
ments of magnetic adatoms on metallic surfaces. 
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APPENDIX A: SIMPLIFIED EXPRESSIONS FOR 
THE FIRST THREE NONLOCAL SURFACE 
GREEN FUNCTIONS 

In this appendix, we provide simplified expressions for 
the nonlocal surface Green functions Goi, Gi%, and G02, 
involving just a single integration. These expressions are 
analogous to Eq. ( p7| ) for Goo- 

Denoting for convenience z = e + if] and Cfe — z ~ 
2tcos(fc), the nonlocal Green functions are expressed as 



with 



dk 

G lm {z)= I G 2D (( k )F lm (z,k) — , 



it' 

a 



Fqi(z, k) — cos(fc) — cos(2fc) 



(Al) 



(A2) 



F n (z, k) — cos(fc) 



cos(2fc) 



2 6i 2 



z 
3t 



z(k 
12< 2 ' 



(A3) 



and 



-F02O, k) = cos(fc) \ [ - 1 



it 



2z 
St 



cos(2fc) 1 2 



12i 2 



cos(3fc) 



a' 



(A4) 



Similar expressions, but with modified Fi m (z,k), apply 
also to Gi m (z) with larger values of \l\ + \m\. 



Refs. |] and ||, which focused on the case where Ri = R s . 
In Ref. [l], Li et al. considered the case of zero direct 
tunneling between the tip and adatom, corresponding to 
q = in Fano's notation. The antiresonance that de- 
velops in G(V) in this case was approximated by the in- 
verted NCA lineshape of the impurity spectral function, 
computed within a degenerate Anderson model that ac- 
counts for the full 4/ degeneracy in Ce. The depth of the 
antiresonance was left as a fitting parameter. 

As evident from Eqs. (|^)-(|lO|) with Ri — R s , the above 
relation between the impurity contribution to the differ- 
ential conductance and the impurity spectral function is 
exact in the limit T — ► 0, provided Goo(e + if]) = —i^Po 
is essentially constant for |e| on the scale of Tk- This 
relation loses accuracy, though, for T > Tk, when the 
convolution with the derivative of the Fermi-Dirac func- 
tion in Eq. (|ic| ) increasingly smears the lineshape of the 
impurity spectral function. 

Contrary to Li et al., who restricted attention to q = 0, 
Madhavan et al. considered the full range in q. To this 
end, Fano's expression for the differential conductance 
was generalized according ton 



G(V) = G, 



(q + eV) 2 
1 + (el/') 2 ' 



eV' = 



eV-ed-Re^jjeV + ir])} 
Im{S^(eV + iri)} 



(Bl) 



(B2) 



where (e + iff) is the full d-electron self-energy, includ- 
ing both the on-site repulsion U and the hybridization to 
the conduction band, Vh- The d self-energy was approx- 
imated in turn by a form corresponding to a Lorentzian 
Abrikosov-Suhl resonance, with a half-width Tr and a 
peak position that was left as a fitting parameter. 

Compar ison with Eqs. (||)-([l0]) for Ri — R s reveals 
that Eqs. ( |Bl[ )-(p^) are correct in the limit T — > 0, pro- 
vided Im { S^(e + irf) } = —T. The latter equality is exact 
for T = and e = 0]l3 and is a reasonable approxima- 
tion for T < Tk and |e| ^ Tk- This approach, however, 
breaks down for T > Tk, both due to the inapplicability 
of the assumed form of X^(e + irf), and because of the 
convolution with the derivative of the Fermi-Dirac func- 
tion in Eq. 
of Gf a (e + iijy. 



(10) which smears the underlying structure 
Hence, similar to the analysis of Li et at, 



this approach is restricted to the low-temperature regime. 
By contrast, Eqs. (p|)-(|l0|) are valid for any temperature 
T and any R s ^ Ri, and are easily amendable [using 
Eqs. (|23|)1 to the case of nonlocal tunneling between the 
STM tip and the substrate conduction electrons. 



APPENDIX B: COMPARISON WITH PREVIOUS 
WORK 

It is instructive to compare the present theory of point 
tunneling between the STM tip and the substrate con- 
duction electrons, Eqs. (j^)— ([Toj) , with the analyses of 



APPENDIX C: SEVERAL MAGNETIC 
ADATOMS 

In this appendix, we generalize our formulation of the 
tunneling current to the case of several magnetic adatoms 
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deposited on top of the metallic surface. Specifically, we 
consider a cluster of m close-by adatoms positioned at 
points Rj (j = 1, • • • , to), each with its own hybridiza- 
tion matrix element, Vhj, and its own tunneling matrix 
element, tj. The different adatoms need not be identi- 
cal, and can generally have different <i-level energies and 
different on-site Coulomb repulsions (denoted by tj and 
Uj, respectively). The Hamiltonian of the system has the 
form H = TLsub + Hu P + Htun, where Hu P is described 
by Eq. (0), and H su b and Tt tU n are given by 



H sub = + E <U E n % + U i 

ka 3 y a 

+ E v m + tUfydj*} , (ci) 

■Htun = t c E{^(^ S )^ + 4Vv(-R S )} 

<7 

+E^E{4^+4^}- (C2) 

3 a 

Here dt creates an atomic d electron with spin a on 

the j-th adatom, and n^ CT = dj^dja is the corresponding 
number operator. All other notations are the same as in 
Eqs. 

Evaluating the tunneling current from the substrate to 
the tip to second order in t c and tj, one obtains an ex- 
pression identical to that of Eq. (0), with the sole mod- 
ification that /0/o-(e) = — ImG/o^e + it]) represents now 
the zero-tunneling spectral function corresponding to 



fa = t c 1p a (R s ) + E *i d i<r< 



(C3) 



Introducing the m x m matrix d Green function 

/OO 
Gf^t'y^-t'Ut, (C4) 
-oo 

G% a (tX) = -e(t-t')({d ia (t),dl(t')}), (cs) 

together with the two "vector" quantities 

Vj (e + irj) = tj + l,\ h! (;, ; l} (e + i V ), (C6) 



u 3 (e + ii]) = tj + t c VhjGft,^(e + irj), 



(C7) 



the retarded / Green function is conveniently expressed 
as 

Gfaie + irj) =t 2 c Gfi a j ia (e + ir)) 

+ Vi(e + iv)Gt j(T (e + iv,) Uj {e + irj), (C8) 

i,3 

which has the compact matrix representation: 

G fa (e + ir))=t 2 c G MsA (e + ir l )+ [^Gju] (e + »»/). 

(C9) 



All information of the adatom cluster and its many- 
body physics is contained within the G f„ Green function 
of Eq. (C8), which replaces that of Eq. (^) in the final 
expression for the differential conductance, Eq. (|Io|). In 
particular, Eq. (j^) is properly recovered in the case of 
just a single magnetic adatom. This permits the analysis 
of complicated multiple-adatom config ura tions in terms 
of the matrix d Green function of Eq. (C4). 
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